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We consider as a clue to the understanding of the pseudogap phase in High T c superconductors the 
metal-insulator crossover in underdoped, non-superconducting cuprates as temperature decreases 
and a similar crossover in superconducting cuprates when a strong magnetic field suppresses su- 
perconductivity. A spin (SU(2)) and charge (U(l)) Chern-Simons gauge field theory, applied to 
the t — J model is developed to describe this striking phenomenon. Two length scales have been 
derived from the theory: The antiferromagnetic correlation length £ « (<5| In <5|) _1/ ' 2 , where S is 
the doping concentration and the thermal de Broglie wave length of the dissipative charge carriers 
At ~ (TS /t)^ 1 / 2 , where T is the temperature, t the hopping integral. At low temperatures £ < At, 
the antiferromagnetic short range order dominates, and the charge carriers become localized showing 
insulating behavior. On the contrary, at high temperatures £ > At, the dissipative motion of charge 
carriers prevails, exhibiting metallic conductivity. Furthermore, the gauge interaction induces bind- 
ing of spinon (spin excitation) and holon (charge carrier) into "electron" resonance. This process 
introduces a new energy scale, the inverse recombination time, which turns out to be essential in 
the interpretation of the out-of-plane resistivity. The major steps in the theoretical derivation of 
these results, particularly the calculation of the current-current correlation function and the Green's 
function for the physical electron, are presented with some detail. The obtained theoretical results 
are systematically compared with the in-plane and out-of-plane resistivity data, the magnetoresis- 
tance data, as well as the nuclear magnetic relaxation data in the pseudogap regime. A very good 
agreement is obtained. 

PACS numbers: 71.10.Hf, 11.15.-q, 71.27.+a,74.25.Fg 

I. INTRODUCTION 

The mechanism of High T c superconductivity, a great intellectual challenge posed by Bcdnorz and Miiller to the 
condensed matter physics community 16 years ago, to large extent remains unresolved. The understanding of the 
"pseudogap", i.e., a suppressed density of states (DOS) near the Fermi level, as one of the most distinctive features 
of cuprates (particularly in the underdoped regime) from the conventional superconductors, may serve as a key to the 
High T c conundrum. ( For a recent review on this issue please consult Ref. Q.) 

A. Pseudogap Phenomenology 

A gap in the spin excitation spectrum was first observed in the Nuclear Magnetic Resonance (NMR) experiments^ 
showing up as DOS reduction in both Knight shift and spin-lattice relaxation rate, below certain crossover temperature 
T*. Later on, a gap in the charge excitations was also found in transport measurements. Unlike optimally-doped 
cuprates where the linear temperature dependence of the in-plane resistivity extends to very low temperatures, in 
underdoped samples the temperature dependence becomes sublinear below some characteristic temperature T* because 
of reduced scattering rate due to gap opening^^ Upon further decrease of doping a striking metal-insulator crossover 
(MIC) was observed^ The existence of pseudogap in the normal state of underdoped cuprates was further confirmed 
by studying the a&-plane optical conductivity^ the specific heat measurements ji and many others. However, the most 
direct evidence of the pseudogap comes from the Angle-Resolved-Photoemission Spectroscopy(ARPES) experiments 
which showed a rf-wave like gap in the normal state in Bi-2212 below the crossover temperature T*A The pseudogap 
structure in the total density of states was also detected in the tunneling experiments on Bi-2212i2 

The current picture of the pseudogap opening appears as follows^ There exists a large Fermi surface (FS) in 
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cupratc superconductors, consistent with predictions of electronic structure studies at high temperatures. As the 
temperature decreases below T* , the pseudogap first opens near (0,7r), it then gradually "eats" up the original FS, 
converting it into "pscudogapped" part, eventually leaving only short disconnected arcs around (7r/2,7r/2). Finally, 
these arcs shrink to nodal points of the gap function, and the pseudogap converges to the superconducting gap, with 
the same rf-wave like symmetry. Meanwhile, the quasiparticle peak which is ill-defined in the normal-pseudogap state, 
becomes well-defined in the superconducting state, like in the typical BCS superconductors. This picture seems to be 
universal for all classes of cuprate superconductors, and is consistent with all available data on pseudogap phenomena. 

B. Different Schools of Thought 

However, the smooth crossover of pseudogap into superconducting gap does not tell us the origin of the pseudogap 
itself. Theoretical models attempting to interpret the pseudogap are very much diversified, including the nearly 
antifcrromagnetic Fermi liquid (NAFF) approachjii the crossover from Bose-Einstein condensation(BEC) to BCS 
scenario^ and many othersi Closely related to our approach, there are two schools of thought. According to one 
school, pseudogap is a "precursor" to superconductivity, i.e., the superconducting state is more fundamental. The 
other school emphasizes the proximity to the Mott insulating state, considering the pseudogap phase as nothing but 
doped Mott insulators. Therefore, the normal state is more fundamental, while the superconduting state is derived 
from this anomalous normal state, as the BCS superconductivity appears when a pairing force is present in the Landau 
Fermi liquid. 

There were several early proposals on the first approach, assuming preformed pairs in the normal stated This 
approach was developed further by different groups. Emery and Kivelson considered phase fluctuations of the su- 
perconducting order parameter which destroy the coherence above T c ^4 Somewhat related appear the nodal liquid 
approach^ and its QED3 variant. Randeria and collaborators attempted to derive several properties of the pseudo- 
gap phase, starting from a low-density, short-correlation length superconducting stated It's fair to say that to describe 
the rich variety of phenomena in the pseudogap phase by superconducting fluctuations alone is a too difficult, if not 
impossible task. 

The other school of thought was pioneered by P.W. Anderson^ The reference compound before doping is a Mott 
insulator characterized by a strong on-site Coulomb repulsion. The electrons are basically localized, forming singlet 
pairs. A finite amount of energy (spin gap) is needed to break this singlet pair. Although similar to superconducting 
singlet pairing, this state itself is very different from the superconducting state. A mean field theory using the 
spin-charge separation concept and "slave boson" technique was developed very early^S To implement the single- 
occupancy constraint coming from the on-site repulsion and to describe the interaction between "slave" particles a 
gauge field theory was developedi 20 i 21 First the U(l) gauge field theory^ then the SU(2) gauge field theory 2 ^ was 
developed to account for various properties of the pseudogap phase by P.A. Lee and his collaborators. The transport 
properties in the pseudogap phase was also considered using the U(l) gauge field models A new version of this 
approach was formulated recently in terms of "spinons" , "chargons" and "visons" with some predictions to be verified 
by experiments^ 

A different, but related approach is to associate the pseudogap phase with a quantum critical point (QCP)^ It was 
proposed that the strong fluctuations near the QCP may be responsible for superconducting pairing and anomalous 
properties in the normal stated What is the order parameter related to this quantum phase transition? There 
were several proposals 2 ^ 2 ^ to be checked by experiments. By extrapolating the curve T*(6), where S is the doping 
concentration to temperatures below T c , a "critical concentration"^ ~ 0.19 was identified, where the superconducting 
condensation energy also reaches a maximum^ The specific heat jump at the superconducting transition changes its 
doping dependence at this point as wellS These authors argue that below the QCP the short range antifcrromagnetic 
correlations dominate, giving their way to superconducting ordering at that point. 

Before resorting to review our own work we would like to mention a recent scanning tunneling microscopy (STM) 
experiment exhibiting explicitly the competition and coexistence of the pseudogap and superconducting phases^ The 
STM pattern seems to be messy: the superconducting regions are separated by pseudogap areas at nanoscale, but 
are kept coherent as Josephson arrays. However, these two types of regions refuse to mix with each other, as Zaanen 
put it pictoriallyj^i like "oil and vinegar in salad dressing" . Moreover, these two regions behave totally differently: 
the typical "resonance" states caused by Ni impurities in the superconducting areas disappear completely in the 
pseudogap phase. This means these two phases, in spite of their apparent similarity, arc of different nature^ To 
elucidate their competition and coexistence is a real challenge to theory. 
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C. Clue to the Problem 

What is the key to the understanding of the pseudogap phase? An "obvious" answer is the pseudogap formation, 
i.e., to explain how a gap is originated and how the states are filled in as the temperature increases. Apparently 
this is not enough, as we know that one scale or a single variable usually cannot give rise to a "show" in physics. 
There must be some competition. Our attention was attracted by the spectacular MIC observed in undcrdopcd, 
non-superconducting cuprates in the absence of magnetic field and a similar phenomenon in superconducting samples 
when an applied strong magnetic field suppresses the superconductivity. 

First of all, this crossover is a rather universal phenomenon. A minimum in resistance (around 50 - 100K) 
and a crossover from metallic conductivity at high temperatures to insulating behavior at low temperatures has 
been observed in heavily underdoped L^-xSr^CuCU (LSCO)^i non-superconducting E^+zS^-yCuOeiaj^ and 
non-superconducting YBa2Cu307-,5(YBCO)^^ and La-doped Bi-2201i2I It has also been observed in electron- 
underdoped Nd 2 -2;Ce x Cu04 (NCCO)2& and Pr 2 _ x Cc x Cu04 (PCCO)i^ Very recently, this issue has been studied 
systematically again on a series of high quality LSCO samples^ and the earlier results have been reconfirmed. Quite 
some people consider this crossover as localization due to inhomogeneities. We disagree with this interpretation. In 
some of these samples, the estimated kpl < 0.1, where kp is the Fermi wave vector and I is the elastic mean-free 
path, i.e., the resistivity is well above the Ioffe-Regel limit. This means localization due to disorder effect is irrelevant 
here. 

More importantly, such a MIC has also been observed in a number of superconducting samples when a strong 
magnetic field suppresses the superconductivity. It was first found on LSCO, where a strong pulsed magnetic field 
up to 60T could suppress superconductivity in samples up to optimal doping^i^ Further, a crossover to insulating 
behavior was found at low temperatures. Similar crossover was also found in La-doped Bi-2201^i4 electron-doped 
Pr2_ 2; Ce :r Cu04^ as well as in Zn-doped YBC04^ Again, the interpretation in terms of localization due to disorder 
does not work here, since the estimated kpt ~ 12 — 25 at the MIC point in some of these systems, and the resistivity 
is well below the Ioffe-Regel limit. Ascribing MIC to proximity to the QCP as caused by charge density wave or 
d-density wave instabilities 2 ^ 2 ^ is not a good explanation, cither. While in LSCO the insulating behavior persists up 
to optimal doping 41 i 42 in La-doped Bi-2201 such a behavior stops at 1/8 doping^ well below the optimal doping, and 
no signatures of any stripe phase showed up. Several authors consider the insulating behavior at low temperatures 
as due to non-Fermi liquid properties^ However, the MIC itself was not addressed in any theoretical considerations 
known to us. 

In our view the MIC is the clue to the understanding of the pseudogap phase. We consider the MIC in underdoped 
cuprates in the absence of magnetic field and MIC in superconducting samples when a strong magnetic field suppresses 
superconductivity the same phenomenon with the same origin: as an outcome of competition between the short range 
antifcrromagnctic order and the dissipative motion of the charge carriers. We have applied the SU(2)xU(l) Chern- 
Simons (the spin-charge) gauge field theory to treat this problemi 4 ^ 4 ^ 4 ^ The formalism itself will be outlined briefly 
in the next Section, whereas here we give some intuitive picture of our main results. 

We start from the Mott insulating state which shows antifcrromagnctic long range order (AF LRO). Upon doping 
beyond certain threshold the AF LRO is destroyed, being replaced by AF short range order (SRO), characterized 
by an AF correlation length £. Since the holes distort the AF background and their average distance ~ <5~ 1//2 , 
intuitively, £ sa (J -1 / 2 , where S is the doping concentration. This has been confirmed by the neutron scattering 
experiments^ Our theoretical treatment proves £ ks (5\ In <5|) — 1 / 2 , providing the first length scale. Put another way, 
the corresponding energy scale is the spin excitation (spinon) gap m s — J(S\ ln^l) 1 / 2 , where J is the AF exchange 
interaction. A competing factor is the dissipative motion of the charge carriers with characteristic energy ~ Tm^, 
where T is the temperature, while mh ~ 5/t is the effective mass of the charge carrier (holon), t the hopping integral. 
The corresponding length scale is the thermal de Broglie wave length At ~ (TS/t)^ 1 ' 2 . We know that in transport 
and many related phenomena when two scales are competing producing a crossover, only the shortest time (and 
corresponding length) scale, or the largest energy scale matters. Thus at low temperatures £ < At, the AF SRO 
dominates, the charge carriers become localized showing insulating behavior. We would like to emphasize this "peculiar 
localization" is mainly due to interaction rather than disorder. On the other hand, at high temperatures £ > At, 
the dissipative motion of charge carriers dominates, exhibiting metallic conductivity. Therefore the competition, or 
fighting of the real part of the "self-energy", the mass gap, and the imaginary part, the dissipation gives rise to this 
spectacular phenomenon: MIC. In this sense, it serves as a clue to the understanding of the pseudogap phase. 

A natural question would be: what are the other consequences of the existence of these two factors? It turns out 
that a number of experimental observations in the pseudogap phase can be explained by it. Roughly speaking, the spin 
(spinons) and charge (holons) excitations in the pseudogap phase behave like "separate particles" in their scattering 
against gauge fluctuations, which renormalizes their properties and dominates the in-plane transport phenomena. 
However, at small enery- momentum scale the gauge field binds spinon and anti-spinon into magnon "resonance", 
similarly spinon and holon are bound into electron "resonance" with non- Fermi-liquid properties. In particular, these 
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" recombincd" particles show up in the out-of-plane transport. 

D. Outline of the Rest of the Paper 

The formulation of the spin-charge gauge field theory was presented earlier^ The calculation of the in-planc 
resistivity and MIC in underdoped cuprates in the absence of magnetic field was briefly reported in Ref. l48[ whereas 
the MIC in superconducting samples when a strong magnetic field suppresses superconductivity was considered in 
Refl49l Due to space limitations the presentation in the these two publications was inevitably too concise. The main 
purpose of this paper is twofold: To present the computations in more detail to outline the essential steps and to 
compare the results of our calculation with experiments in the pseudogap phase in a more systematic way. Most 
results have not been published before. 

The spin-charge gauge formalism is outlined in Section II, whereas the calculation of the in-plane and out-of- 
plane resistivity is described in Section III. The computation of the "spinon" current-current correlation function 
which is the key ingredient in studying many physical quantities is presented in Section IV. The calculation of the 
Green's function for the physical electron needed to compare with the ARPES data and FS, as well as the c-axis 
resistivity calculation, is described in Section V. In our approach the spin and charge are not fully separated, as in 
one-dimensional interacting systems. They arc not confined, either. Instead, they form a bound state due to the 
transverse gauge field. A new energy scale, the inverse recombination time for the physical electron, appears in the 
binding process which shows up in the out-of-plane resistivity. Section VI is devoted to comparison of theory with 
experiment. We first start with the in-plane resistivity (A), whose normalized value shows universal behavior (B), 
then continue with out-of-plane resistivity (C). Finally, wc consider other observables, including the magnctorcsistancc 
(D) and the spin-lattice relaxation rate (E). The paper ends with several concluding remarks (Section VII). 

II. THE SPIN-CHARGE GAUGE FORMALISM 

In this Section we review the formalism involved in the derivation of the low-energy effective action for the t-J 
model in a range of parameters which should provide an adequate theoretical description of the "pseudogap phase" 
of high-T c cuprates. This effective action will be the starting point for the computation of physical observables to be 
compared with the experimental data. For further details on the derivation, see Ref. l47l 

A. Chern-Simons Representation of the t — J Model 

Our theoretical treatment of the t-J model is based on the following theoremt^^i 

If we couple the fcrmions of the t-J model to a U(l) gauge field, B^, gauging the global charge symmetry, and to 
an SU(2) gauge field, V^, gauging the global spin symmetry of the model, and we assume that the dynamics of the 
gauge fields is described by the Chern-Simons actions: 

Sc.sXB) =~J (fxe^B^Bp, 

Sc.sXV) = ±-J d 3 xTi-e^[V^V p + ^V^Vp], (1) 

where & vp is the fully anti-symmetric tensor, then the spin-charge (or 5/7(2) x gauged model so obtained is 

exactly equivalent to the original t-J model. 

Let us give some ideas of the proof of the above theorem for the partition function. We expand the partition 
function of the gauged model in the first-quantized formalism in terms of the worldlines of fermions. After integrating 
out the gauge fields, the effect of the coupling to B^iV^) is only to give a factor e~' l ^(e 1 ^) for any single exchange of 
the fermion worldlines, so the two effects cancel each other exactly. 

To the fermion field of the gauged model, denoted by Xa (ct spin index), we apply a formal spin-charge decomposition: 
X a ~ Hz a , where H denotes a spinless fermion (holon) field and z a a spin i hard-core boson (spinon) field satisfying 
the constraint z* a z a = 1. This constraint eliminates double occupation, as required in the t-J model. 

The above spin-charge decomposition introduces a further U(l) gauge symmetry which will be called h/s (from 
holon/spinon): 

z a {x) -► e lA{x) z a (x), H(x) -> H(x)e- lA{x \ (2) 
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with A a real gauge parameter, with which a self-generated gauge field is associated, analogous to the one appearing 
in the slave boson and slave fermion approaches 

We remark that, in view of this residual h/s gauge interaction, the spin-charge separation performed above is a 
priori purely formal and only the dynamics of the coupled system determines if it has a physical substance. As an 
example, a confining dynamics would completely destroy the physical spin-charge separation. 

In a Mean-Field- Approximation (MFA)^£ to the spin-charge gauged t-J model, in a region of parameters which 
should correspond to the "pscudogap phase" of high T c cuprates, the role of the three gauge fields is the following: 

- Bmfa carries a flux tt per plaquette, converting via Hofstadter mechanism the spinlcss holon H into a Dirac 
fermion, with linear dispersion and a pseudospin structure related to the two Neel sublattices. The fermion system 
exhibits a "small" FS with ep ~ tS, 6 being the doping concentration, centered around the points (±7r/2, ±7r/2) in 
the Brillouin zone. 

- Vmfa dresses the holons by spin vortices of opposite chirality in the two Neel su blattices. The spinons in the 
presence of this gas of "slowly moving" dressed holons acquire a mass m s ~ \/6\hx5\ yielding SRAFO. This is due 
to a coupling at large scales of the form iyMFA z a z a)- Self-consistency of this treatment relies on the inequality 
cf ~ tS -C e s ~ J \J S\ In S\ for small 8. The derived doping dependence of the AF correlation length is consistent with 
the neutron data^ 

- The self-generated "photon" field couples the Fermi liquid of holons to the gapped spinons, described by a 
massive (CP 1 ) non-linear sigma (NLcr) model. 

B. Magnons and Electrons as "Composite" Particles 

A low-energy effective action for A is obtained by integrating out spinons and holons, in a path-integral formulation. 
We make the assumption that the scaling limit (large distance, long time) can be taken separately for the two 
subsystems^- Then, using the techniques of Ref. 153 . one can prove that in this scaling limit the action is quadratic in 
A. This conclusion follows from a derivative expansion for the spinon action, due to the presence of a mass scale, m s , 
and from a tomographic decomposition along rays perpendicular to the FS of holons, using the quadratic dependence 
on A of the scaling action for a single ray (Schwinger action). This means that all results derived from the renormalized 
quadratic action in A are valid beyond the standard pcrturbative treatment. Thus the skeptism towards the gauge 
field approach based on the worry that the coupling to the gauge field is strong, while the treatment is perturbative, 
is not well justified. 

The term obtained from spinon integration is a 2D Maxwell- like action (as in quantum electrodynamics), because 
the spinons are massive and the spinon action is parity- invariant. The transverse component would then generate a 
logarithmic confining potential between spinons and anti-spinons. The longitudinal part is gapped due to the plasmon 
effect at finite T . This means if only spinons were coupled to the gauge field, the renormalized gauge field would 
have confining dynamics. However, there are also fermions (holons) coupled to the gauge field, as well. The term 
obtained from holon integration, due to the presence of a finite FS, exhibits a Rcizer singularity^ More precisely, the 
transverse component A T of the gauge field turns out to be, for u>, \q\,u/\q\ ~ 0, of the form 

(A T A T )(u } ,q) ^ (- X \q\ 2 + (3) 

where \ i s the diamagnetic susceptibility and k the Landau damping. 

This behavior dominates over the Maxwcllian term at large scales, destroying confinement. Nevertheless, as we shall 
see, the attraction generated by A T in spinon-antispinon and spinon-holon pairs will be sufficient to produce resonances 
with the quantum numbers of the magnon and electron, respectively. Therefore a true spin-charge separation is 
not realized in our approach. An "intermediate" situation in-between confinement and full separation, namely the 
"composite" nature of magnons and electrons is at the root of our interpretation of the anomalous behavior of physical 
quantities. 

A key and novel feature of our approach is the mass of the spinon with a specific doping dependence described 
above. This feature is not shared by the other SU(2) x U(l) gauge field theory2ii where the gauged SU(2) symmetry 
is an enlargement of the particle-hole symmetry at half-filling with switched statistics of holon and spinon w.r.t. ours. 
The mass of the spinons in our approach and its competition with dissipation of the gauge field due to coupling with 
holons have far-reaching consequences, and it turns out to be responsible, in our scheme, for phenomena like the MIC, 
the low-T positive transverse in-plane magnetoresistance, the peak in the DC conductivity and the Cu spin-lattice 
relaxation rate, hence for many experimental signatures of the "pseudogap phase" . 
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C. Motivation for the Gauging Group Choice 



We end this Section with some comments on our choice of the gauging group in the theorem stated at the beginning 
of this Section. This theorem is a special case of the following more general 
ThcoremS&( Chcrn- Simons representations of the t-J model): 

Let G be a subgroup of the global symmetry group of the 2D t-J model; consider the G-gauged t-J model obtained 
by replacing the fermion field, c, of the t-J model in the action, S t j(c), with a new field, x, minimally coupled to a 
gauge field W, with gauge group G. Denote the action of the G-gauged model by Stj(x, W). Define the Chcrn-Simons 
action for W by: 



Then, for a suitable choice of a real constant kc and of the statisics of x, fermionic or bosonic depending on ha, the 
model with action Stj(Xi W) + kGS c . s .(W) is exactly equivalent to the original t-J model. 

The two basic features of the 2D t-J model needed for the proof of the theorem are its dimensionality, necessary to 
apply the Chern-Simons theory, and the Gutzwiller projection forbidding double occupation, needed to have at most 
pointlike intersection between \ worldlines. In fact, (see the comment after the theorem in Section IIA) only with 
this property we can associate well defined phase factors to interchange of the worldlines. 

Each of the Chern-Simons representations allowed by the theorem can be taken as a starting point for a Mean 
Field Theory In particular the slave boson and slave fermion approaches can be derived by choosing G = U(l) and 
&j/(i)=+l and -1, respectively^ Our choice of the SU{2) x U(l) spin-charge gauging is motivated by the analysis 
of the one-dimensional (ID) t-J model in the limit of small J/tmL The model is exactly solvable by Bethe Ansatz 
and Conformal Field Theory techniques, and one knows the critical exponents of its correlation functions. They 
are basically derived by decomposing the fermion field of the model into spinon and holon. The spin-charge gauging 
corresponds, using a dimensional reduction from 2D to ID, to a semionic nature of the two excitations, i.e. an exchange 
of spinon or holon fields yields phase factors e ±l7r / 2 , intermediate between the fermionic e 171 and the bosonic 1. This 
is exactly a property needed to reproduce, in a sort of MFA, the known critical exponents. It may be worthwhile 
to compare the role of the £7(1), SU(2) and h/s gauge fields in ID and 2D. In ID a gauge field has no transverse 
(physical) components, while in 2D it does have one. The disappearance of this degree of freedom in ID w.r.t. 2D 
induces the following effects: 

- Bj IFA =0, hence there is no Hofstatdcr mechanism and the holon has a quadratic dispersion; 

- ^^^=0, hence there is no spinon mass generation; 

- A^=0, hence spinons and holons are decoupled and this yields a true spin-charge separation in ID. 



In this Section we highlight the distinctive features of the experimental data on in-plane and out-of-plane resistivity 
in underdoped cuprates, sketch the theoretical scheme for computations and outline the qualitative understanding of 
the resistivity behavior in our gauge field approach. 



One of the first striking experimental findings on high T c cuprates was the anomalous behavior of in-plane resistivity, 
which in optimally doped samples appear linear in T. In underdoped samples, it deviates from the linear dependence 
at low temperatures, but the standard metallic behavior ~ T 2 derived from the Fermi liquid theory is not observed. 
Instead, there are the following two distinctive features: 

a) In many strongly underdoped samples there exists a minimum in the resistivity, around ~ 50 — 100A', corre- 
sponding to a MIC, as we outlined in Section IC. A similar crossover is also observed in superconducting samples, if 
superconductivity is suppressed by applying a strong magnetic field. 

b) Another characteristic feature of in-plane resistivity which appear quite universal in underdoped samples is an 
inflection point, i.e. a maximum of dp/ dT at T* ~ 100 — 300A; this maximum disappears for higher dopings. At 
even higher temperatures the resistivity exhibits a linear in T behavior approached from below. In the literature T* 
is also defined by some authors as the temperature where the resistivity deviates from the linear dependence, as we 
mentioned in Sec. IA. That value is higher than the inflection point. 




III. RESISTIVITY AND SPIN-CHARGE DECOMPOSITION 



A. In-Plane Resistivity 
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To calculate the in-plane resistivity we use the IofFc-Larkin rule^Si a somewhat counterintuitive but a typical feature 
of the gauge approach, stating that the physical resistivity, p, is a sum of the resistivity due to spinons, p s and the 
resistivity due to holons, p^: 

p = p s +Ph- (4) 

The derivation of this addition rule is based on the following consideration: If wc couple the electron to an external 
electromagnetic (cm.) field, A e . m , it turns out that we can attribute an arbitrary e.m. charge e with < e < 1, to 
the spinon and a charge 1 — e to the holon, because, in the path integral formalism, e can always be eliminated by 
the change of variable A — > A + eA e . m . As a consequence, neglecting "photon" drag, the renormalized e.m. current 
polarization bubble, IT e m , obeys the rule: 

(iwr 1 = (u s )- 1 + (Uh)- 1 . (5) 

From (J2J) and Kubo formula, one can derive Ioffe-Larkin rule, provided both conductivities a s and ah are non- vanishing. 
This will be self-consistently verified a posteriori except for very low temperatures. A crucial assumption here is the 
quadratic dependence of the effective action in A, which, as we pointed out in Sec. IIA, is valid beyond the standard 
perturbation expansion in the scaling limit. So does the Ioffe-Larkin formula. 

Denoting by j s the spinon current, the spinon resistivity is calculated from the fully renormalized current polariza- 
tion bubble: 

(f(x)f(y))=IL s (x-y) (6) 

via the Kubo formula: 

(Ps)' 1 = <T S = -(w- 1 Imnf(w,g = 0))| w ^ 

✓>oo 

= 2 / dx°x°n s (x o ,q = 0). (7) 
Jo 

For holons we have similar equations replacing the index s by h (e.g. j denotes the holon current). In eq.© the 
expectation value is taken by integration over A, the spinon field, z a of the continuum NLer model and the Dirac 
holon field ip. The last equality in JJJ is obtained via Lehmann representation and the superscript R denotes the 
retarded propagator. 

As pointed out first by Anderson^ the in-plane resistivity should be interpreted in terms of spin-charge separation. 
In a gauge approach if the scattering time of spinons or holons by gauge fluctuations is shorter than the lifetime of 
the electron (as in our case), then this time scale will dominate the in-plane resistivity and it might exhibit a different 
temperature dependence than the electron lifetime. 

It turns out that the peculiar feature like MIC is mainly due to the spinon contribution. As we mentioned in 
Sec.IC, MIC is caused by the competition between the spinon mass term with the gauge field dissipation. The spinon 
contribution to resistivity, proportional to the spinon scattering time against the gauge field, turns out to be ~ T _1 
at low temperatures when the spinon mass effect dominates, while it is ~ T 1 / 4 at high temperatures when the gauge 
field dissipation overwhelms. 

In our theoretical framework we identify the inflection point as the signal of a crossover to a different "phase" , the 
so-called strange-metal phase, characterized by T-linear resistivity, which will be addressed in a separate paper^ 
Moreover, as will be shown in Sec. VIB, the normalized resistivity is a universal curve, if its value at the MIC point 
Tmic and the inflection point T* are used as references. At very low temperature many samples exhibit a second 
inflection point below which the resistivity appears approximately logarithmic in T)ii again we interpret this as a 
crossover to a different "phase" . 

B. Out-of-Plane Resistivity 

The out-of-plane resistivity exhibits a completely different T dependence. At low temperatures in the "pseudogap 
phase" p c is insulating, behaving like T~ x with a coefficient essentially independent of the material^ At higher 
temperatures p c typically develops a rounded kneei^iSi As emphasized by Anderson^ the coexistence at the same 
temperature of a metallic in-plane and an insulating out-of-plane resistivity is hard to reconcile within a Fermi liquid 
theory, whereas it might have a natural explanation in the framework of spin-charge "separation" . In such scheme, 
in fact, a spinon- holon decomposition of the electron holds only in the CuC>2 layer and spinons and holons should 



8 



recombine into electrons to hop between layers and contribute to p c . The out-of-planc resistivity is then determined 
by the time scale of electron recombination. In our approach we have a way to implement this general ideas proposed 
by Anderson in the "pscudogap phase" . 

To calculate p c we use the approach proposed by Kumar and Jayannavar (K-J)§i which is motivated by the 
experimental observation that the c-axis transport is essentially incoherent, i.e. there is no band-like motion orthogonal 
to the CuC>2 planes. One can then consider a system of two layers weakly coupled by an effective tunnelling matrix 
element — t c , taking into account an averaged momentum dependence of the hopping parameter (vanishing for diagonal 
momenta). One can write the 2D retarded Green function of the electron (holon-spinon) resonance for small to and 
momentum kp on the FS as 

G R (w,k F ) ~ —j^, (8) 
uj + II 

where Z is the wave function renormalization and T the scattering rate. Taking into account a virtual hopping between 
two layers induces a shift of the real part of the denominator of © from ui to u> ± Zt c . Let us denote by G± the 
corresponding Green's functions. The out-of-plane conductivity in the incoherence regime can be written through the 
Kubo formula as 



E 



^2tle 2 ^{oj)A + (k,oj)A^k,oj), (9) 



where A± = — ^Im G± are the spectral functions and the Fermi distribution function. Inserting (|SJ| in @ after 
standard manipulations one obtains 

Po--.(^ + tLo], (10) 



v \Y tlZ* , 

where v is the density of states at the FS. One can already anticipate that the first term causes the insulating behaviour 
and, being independent of i C! it is essentially independent of the material, as experimentally observed. Via eq. fTDjl 
we have related the behavior of p c to the computation of F and Z, thus to the low-energy behavior of electron Green 
function. This propagator in turn can be expressed at large scales in terms of holon and spinon fields, z a , ip a and be 
extracted from a linear combination of terms 

(iW^WiW^W). (n) 

where a denotes the pscudospin structure of the Dirac holons, a the component-spin index of the spinons, and the 
propagator is calculated using the low-energy spinon and holon effective actions. 

We shall see that the derived lifetime of the electron resonance is ~ y-i j 1 -!/ 2 a ^ i ow anc j high temperatures, 
respectively, and therefore it cannot explain the temperature behavior of in-plane resistivity, in particular the MIC, 
but it indeed sets the scale of the out-of-plane resistivity. We shall also see that the theoretical curve indeed has a 
rounded knee, corresponding to the crossover between the high- and low-temperature regimes, in full consistency with 
experiment. 

The following two Sections (IV and V) are more technical. Those who are mainly interested in the qualitative 
aspect of the gauge field approach can skip them and move directly to Section VI. 

IV. THE SPINON CURRENT-CURRENT CORRELATION FUNCTION 

In this Section we outline the computation of the spinon current polarization bubble, i.e., the current-current 
correlation function II s (w, q). at small u> and q. This computation was briefly sketched in Ref. and it is needed 
to derive the in-plane resistivity, as explained above. We will provide more technical detail here for those who would 
like to follow the actual calculation. 

A. Feynman-Schwinger-Fradkin Representation 

We start by writing explicitly the spinon NLer model effective action 



S = I d 3 x- 
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v~ 2 | (d - iAo) z a \ + \(di - iAi) Za\ + m 2 s z* a z a , (12) 
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where g ~ v s ~ Ja is the spinon velocity, with a the lattice spacing, and mj 1 ~ a/ ( | (5 In 5| ) the spinon 

correlation length. After a suitable riscaling of variables, the spinon propagator can be recast in the Schwinger 
representation: 

/>oo 

G a (x,y\A)=igv s dse^^+^fay). (13) 



where x — (v s X , x), A — (vsAq^A) and Aa denotes the 3D covariant Dalambcrtian (or relativistic Laplacian). The 
propagator has been considered in the zero temperature formalism, an approximation justified by the mass gap of the 
spinon. Roughly speaking, it is valid provided T -C Jam s ~ J(\S In (5Q 1 / 2 . The kernel appearing in H13f) has the formal 
structure of an evolution kernel for a 3D Hamiltonian H = — A ,4 + ra 2 s and time parameter s. It can thus be expanded 
in terms of Feynman paths starting from y at "time" and reaching x at "time" s. It is convenient to parametrize 
these paths through their 3- velocity, /i = 0,1,2, using a Feynman-Schwinger-Fradkin (FSF) representation (see 
e.g. Ref. |H: 

G a (x,0\A) = igv s dse- ism > f V<f> [ d 3 pe ip(/ ° *(*)«*-*) e */o* «a[i^(t)+^A(»+/ V(f )««')]. (14) 



Here the p-integration enforces the constraint on the initial and final points of the paths and we use a short-hand 
notation for the 3D scalar product: e.g. 

p-x=p ft x fi . (15) 

For a better understanding of the formula Q14|l , notice that formally setting </) M (t) = dx d ^ the last exponential is i 
times the Lagrangian of a 3D particle coupled to the e.m. potential A^, corresponding to the previous Hamiltonian 
-ff, as one expects in a path-integral formulation. Since under a h/s gauge transformation A(x), the spinon field z a [x) 
changes by the phase factor e lAt - x \ it follows that 

G a (x,0\A„ + d ll A) = e i ^- A ^G ct (x,0\A ll ). (16) 

The gauge dependence of the Green function is already captured by the so-called "Gor'kov approximation" 

G a (x,0\A) = e^fo A ^G a (x,0), (17) 

where J Q X denote integration along a straight line from to a; and G a (x,0) is the free propagator (in the absence of 
gauge field). The expression (|14|) is useful to go beyond Gor'kov approximation by means of the identity^ 

A„(x+ f </)(t')dt')<t> fl (t)dt= f A^dx"- 



dXX f dt 
Jo Jo 



dt'^(t)r(t')F^ fx + X cj>(t")dtA , (18) 



where F^„ = d^A v — d^A^ is the gauge field strength. Pictorially, this representation is illustrated in Fig. ^ The 
second term in (|18|l denoted by S(-P), gives the correction to Gor'kov approximation and it is gauge invariant, as it 
depends only on F^. Shifting 4>^{t) by 2p^ one can rewrite 

G a (x,0\A) = e l ^ A ^G a (x,0\F), (19) 



/ e -i P x-i(p 2 +m*) 



(20) 
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B. Gauge Field Strength F Correlation Function 

Now we turn to the polarization operator Il s . Expressing it in terms of spinon propagators we find: 

U s (x, y) = (D A(x) G(x, y\A)D\ {y) G{y, x\A)) A = ((fy - l -f F^)G{x, y\F)(d» - % -f F^dx p )G(x, y\ - F)) A , 

(21) 

where (■} A denotes the integration over A with the effective action in the scaling limit and D A the covariant derivative. 
Notice that the two non-gauge invariant Gor'kov terms of the two spinon propagators cancel against each other so 
that the result is explicitly gauge-invariant and it depends only on F. We use now the quadratic structure of the 
scaling action S(A) (see Sect. II B) to integrate out the gauge field. The explicit expression for S(A) in the Coulomb 
gauge is : 

S(A) = ^ J dx°(fxA^W v A v , 

with non-vanishing polarization components in the limit of small w, q and ui/\q\, where one finds the leading, Reizer 
singularity, given by^ 

Ui j (q,u;) = (S ij -^)l-iK^+x\tf], i,j = l,2, (22) 



U 00 (q,uj) = v + lo p . (23) 

In H22fl and 123|) x — Xs + Xh, where Xs(h) i s the spinon (holon) diamagnetic susceptibility, k the Landau damping, 
v the density of states at the FS of holons and ui p the plasmon gap. For free holons 

^ H 127TTO/j 67T<5 ^ 

and for free spinons 

Xs ~ m-\ (25) 

Hence, in this presumably reasonable approximation, for low doping Xh » Xs- Due to the dependence on the field 
strength F in (|21|l only correlator of electric (-Fbi) and magnetic (i^j) field can appear in the computation. Since the 
Ao-propagator is short-ranged whereas the A T -propagator is long-ranged, the "electric" field contribution at large 
scales is negligible w.r.t. the "magnetic", and in first approximation we neglect it. However, it might be useful to 
keep in mind that doing this we neglect a short-range attraction between spinon and antispinon (or holon). Due to 
the gapless nature of A T we consider the effect of finite temperature using the thermal propagator: 



duj f dk \k\ 2 e - t ^"+ lk - 



(F ij{ x)F rs (0)) = {S ir S js - 5 is S jr ) J - J — >> ^ 3 coth(-). (26) 

This does not contradict our earlier approximation in considering the zero-temperature spinon propagator in view 
of the finite spinon mass gap. The leading order in T correction enters via the thermal gauge-field propagator in our 
scheme. Since the energy scale for field fluctuations is set by T, in l|26() the integration over frequency is cut-off at 
uj <T which in turn implies \k\ < (^ 1 ) 1 ^ 3 . 

In the limit T£° < 1 an approximate evaluation of the above integral gives 

-i—Qle-^, (27) 

1 /3 

where Q a = (^^^j is a momentum cutoff and Qq 1 can be identified as the length scale of gauge fluctuations, 

analogous to the anomalous skin depth. [^4 posteriori the upper limit for T in the inequality above turns out to be 
reasonable because the typical time scale is ~ 1 ~ T -1 / 3 < T _1 at low T] It turns out that this scale is triggering 
also the size of the spinon-antispinon magnon resonance. It also follows that for m s Qq in the expectation value 
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(|21|l the derivative term dominates over the i^-terms at large scales, so that to evaluate II s the leading term is obtained 
by computing 

(G(x,0\F)G(x,0\-F)) A (28) 

and then taking the spatial derivatives. Notice that (|28|l coincides with the propagator (fl(x) ■ O(O)), where = z*az 
is a "magnon" field. We denote by 4>^ and </>2 the velocity fields relative to the FSF representation of the two Green 
functions in 1(28(1 . Integrating over A the product of the two FSF expansions one obtains an effective action, I((f>i, 4>2 ), 
in the velocity fields, which is quartic neglecting the ^-dependence in F. This approximation can be self-consistently 
justified a posteriori, because (see eq. <(33ll ): 

p ~ m s > - [' ds'M) ~ s" 1 / 2 ~ (x^/m,)- 1 ' 2 . (29) 
s Jo 



C. Eikonal and Saddle Point Approximation 



The (^-integration is then performed using the eikonal approximation: 

MPH^J# + {J*W<^>*#><*+»«'-, (30, 

where ,a 2 denotes the average w.r.t. the gaussian measure appearing in the l.h.s. of (|30|l . This can be justified if / is 
small, since / ~ T for T small. Within this approximation the contribution of the two Green functions factorizes. This 
factorization in a diagrammatic language means that after the cancellation of self-energy and vertex rcnormalization 
implicitly involved in the cancellation among Gor'kov terms, the remaining leading effect of A-fluctuations is a self- 
energy rcnormalization of the gauge-invariant spinon propagator. At this stage the correlator (f2(x) • O(0)) can be 
written as 



d 3 p / dse { ~ l{p2+m2s ~x na))s+lpx ~x Q ° s29{a)} 



o 



2 



(31) 



where a = \p\sQo, f and g are functions which summarize the effect of gauge fluctuations and their length scale 
variation is in fact ~ Qq . Explicit integral representations of / and g are: 



/(a) = a 2 [ dXX f dXX f dv ^-^(a-a) 2 ^ 
Jo Jo Jo 

g(a) = f dXX f dXX f dwe'^^-V 2 . (32) 
Jo Jo Jo 

2 > 



Finally we evaluate the p and s integrals by saddle point approximation, obtaining for m s > T/\- 



1 (x°) 2 -f2 

p~x/2s s~- \> (33) 
V s ~ x^°^ 



and the magnon (l propagator in x-spacc becomes 



(0W ' n(Q)) ~ Wr^ e ( 34 ) 

To apply the Kubo formula one needs to perform the Fourier transform at q = of (j s (x)-j s (0)) ~ (d^fl(x)- 8^^1(0)). 
We consider the region x° >> \x\ and evaluate the | x | -integration via saddle point. Using the form of / and g one 
finds that the exponent of 1)34(1 at large a; exhibits a complex saddle point at a scale (a? ) ~ (x ) 1 / 2 , thus verifying 
the above assumed inequality with a behavior of a standard diffusion, and with argument 7r/4. [In the more precise 
numerical evaluation we neglect small scale fluctuations splitting the above saddle-point into a set of isolated saddle 
points.] A numerical extrapolation in the region of small x° yields an approximate cc°-dependence of the form: 

\x\(x°) ~ e™/ 4 x c (x°), x c (x°) = (C 2 Q^ 2 + C'\x \/m s ) 1 / 2 (35) 
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with C,C finite positive constant (C ~ 0.5), thus approaching a finite value as x° — > 0. Setting a(x°) = Qq\x\(x°) 
and 

= -Hlml ~ IfWx ))* ~ ^Q 2 o9 (a(x )) { ff (36) 
V X 4x m*- ±f(a(x )) 

we have: 

(j^)(<f=0,x°) = J d 2 x(j^x° 7 x)^(0 7 6)) 

~ ^ r ^ ~ 1/2 e -^°) (37) 

((x ) 2 -a;2)3 ^^(x ) 2 / ' { ' 

Since / is smooth on the scale ~ Qq 1 , assuming for x° the same scale the dominance of the saddle point requires 
a lower bound for the temperature, which combined with previous upper bound yields a range of validity given by 

m 2 s >-> m s Q . (38) 
X 

In physical units, this gives a range of temperatures between a few tens and a few hundreds of Kelvin. The real part 
of the exponential in l|37|) is monotonically decreasing in x°, therefore we evaluate the x° integral appearing in Kubo 
formula J7J by principal part evaluation. Since our approach is valid only at large scales, we introduce an UV cutoff 
in the integral at AQq 1 and evaluate the integration assuming A large. Then we make the conjecture that for small uj 
the physics is dominated by large scales and the small-scale contribution can be taken into account by removing the 
UV cutoff after a multiplicative scale renormalization. The result of this approximation is 

er s = 2 Im / dx°x°n{x°, q = 0) ~ Im f Zj I , (39) 

Jo \y/m 2 s-$f(Ce™/*)J 

where Zj = QqZq, Z n = (m 2 . - ic|) 1 / 4 (x/T/"(Ce"/ 4 )) 1 / 2 Qj /2 and numerically one finds /(Ce i7r / 4 ) - 0.2 + «3.3 and 

/"(Ce l7T / 4 ) real. For simplicity we set Im /(Ce J7r / 4 ) = c and we still denote by m 2 the quantity m 2 — Re /(Ce"/ 4 )T/x 
which in the range of temperature we are interested is in fact almost equal to m 2 s . 



D. Spinon-Antispinon "Resonance" (Magnon) Propagator 

For a better understanding of the above equation, we notice that the retarded magnon correlator at positive u> is 
given in the same approximations by 

(fi • 6>( W ,$) f a JodflCQo 1 ^ 4 ) (40) 

L) — 2. / m 2 — ic— 

where Jo is the Besscl function. 

The above formula explains the physics underlying ea. (|39|l : the gauge fluctuations couple the spinon-antispinon 
pair into a resonance with mass gap 



T 

mn = 2ReWm 2 — ic— , (41) 



and inverse life-time 



r o 1 = 2Imy m 2 - ic^-. (42) 

Zj and Zq can be interpreted as T-dependent wave-function renormalization factors which modify the temperature 
dependence of the (physical) correlation functions. Finally, we expect that the plausible effect of the neglected residual 
short range attraction is a further renormalization of the mass gap, but we believe that this does not introduce 
significative changes. 
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V. THE ELECTRON GREEN'S FUNCTION 



In this Section we evaluate the continuum limit of the electron Green's function within our approach, extracting, in 
particular, the wave function renormalization constant, Z, and the inverse lifetime, T, needed to compute p c in K-J's 
approach. 



A. Holon Effective Action 



In order to have a more systematic derivation, it is worthwhile to start by writing the hopping Hamiltonian for 
holons, Hhopp, neglecting at first the coupling to the h/s gauge field. 

Restricting the holon field H to the two Neel sublattices, labelled by A, to which the origin belongs, and B, we 
have in momentum space 

- ? M>^>) ( -^il - ~ 2%l o + ) ( £<f ) 

where 

7± = cos^a) ± cos^a), (44) 

a being the lattice spacing and the sum over k running in the reduced Brillouin zone. The eigenvalues of Hhopp are 
given by: 



e±(fc) = ±2fWcos(A; x a) 2 + cos(fc y a) 2 , (45) 

hence they describe double cones with vertices at ^ ± f 7 ±f hi the Brillouin zone. 

Since the chemical potential for the holon system is positive, fi ~ 2tS, only the e+ band of the double cones exhibits 
a FS. For each of these double cones one can identify a two- component, continuum, Dirac field ip a , a =|, J, describing 
the low energy physics of the system. The continuum effective action for holons ip a coupled to the h/s field can 
be cast in the form 

S h (ij>, A) = d 3 x^ Q (d -n- iA ) + v F j t (d t - iAi)]ip, (46) 



where ip = tp'^o, 7o = Pz, 7i = {cy,(T x ) and vf = 2ta being the Fermi velocity. The relation between tp, z and the 
original electron field c a in the two sublattices is found to be given, e.g. near the ^) double cone, by 

<^(*)4») ~ e i; (^^)-((^( a; )^(O)-0 T (^ T (O))z a ( 2 :)z;(O)) (47) 

(cf (aO<£ A (0)) - e<^^ s ((e^^(x)M0) + e^^x)^)) z a (x)z*(0)). 

Analogous relations hold near the other three double cones. Note that in the z correlator, the contribution of the 
spin flips of the "optimal spinon configuration" of Ref. X£\ must be taken into account. 



B. Tomographic Decomposition 

In the previous Section we evaluated the effect of gauge fluctuations on the z correlator at large scales, using 
the FSF path-integral representation. An analogous representation is hard to use for the ip a correlator because of 
the finite density of holons. This representation would in fact contain a series of alternating sign contributions, 
corresponding to an arbitrary number of closed fermion wordlines, descri bing the contributions of the particles in the 
finite-density ground state, besides the path from to x (see e.g. Ref. l5lT) . To overcome this difficulty, we apply 
a dimensional reduction by means of the tomographic decomposition introduced by Luther-HaldaneS To treat the 
low-energy degrees of freedom we choose a slice of thickness A = fcj?/A, with A ^> 1, in momentum space around the 
FS of tp, as shown in Fig. |3 

To simplify the description, we assume a circular FS, an approximation reasonable for low S (the method applies 
nevertheless to the general case by considering a Fermi momentum varying along the FS). We decompose the slice 
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in approximately square sectors; each sector corresponds to a quasi-particle field in the sense of Gallavotti-Shankar 
renormalizationSi (see also Ref. l65j) . Each sector is characterized by a unit vector ft{0), pointing from the center of the 
FS to the centre of the box. labelled by the angle 6 between this direction and the k x axis. The original momentum 
k inside a given sector is written as 



k = k F n(9) 



(48) 



where q spans the box, therefore \q ■ n(6)\, \qA n{9)\ < A. Due to the Dirac structure of ip, to apply the tomographic 
decomposition to the holon propagator, we first decompose the free ip correlator as 



7o(£:o + k F ) + 7^ - iesgn(\k\ ~ k F ) 
1 



-ikx 



d 3 k 
(2tt)3 

d 3 k 

(2tt) 3U k + k F - \k\ + iesgn(\k\ - k F ) 

7o(fc + k F ) - 7 M fc M 
k + k F + \k\ - iesgn(|A:| - k F ) j a[j 



(49) 



In the scaling limit the matrix in square brackets does not have a pole and, for momenta in a box labelled by n{&), 
it approaches 



A(0) 



70 — 7 ■ n(6) 



(50) 



In Ref. [54] it has been shown that the tomographic decomposition is valid at large distances even in the presence 
of a minimal coupling to a "photon" field. Applying the tomographic decomposition to the holon propagator in the 
presence of an external h/s gauge field A, in the scaling limit, using ^49 [l . 1)50 (1 . we derive 



l 



where 



go - H 6i + «esgn(q • n(6»)) 



He = A + n{6) ■ (q- A) + -J- Uq- A) A n{9)] \ 

lk F L J 



(51) 



(52) 



and J A denotes integration over a square box of size A. To the Fourier transform of the term in square bracket of 151|) 
one can apply the FSF path representation (see also Ref. I^for a treatment of linear dispersion). Using manipulations 
analogous to those performed in the previous Section, one can rewrite l|51|l as 



k. 



A 



I deA{6)e- ik ^ s e i ^^ / dqo { 



2ir 



du 



d 2 q 



du 



d 2 q 
{2-Kf 



e(-?n)] 



Jq° (x° -u)-iq(x+n(0)u) 2ft: 



X{ / T>(p±e i ^o ^^J- e -« Jo 1 dTT Io du ' So du"<p^{u')<p ,y (u")F^ v (T f a u v {u"')du"')^ 



(53) 

where we use the short notation qu = q ■ n(9), q± = q A n(0) and ip is the velocity field of components ^(t) = 
(l,l,(pj_(t)). Note that in Q we have replaced the original discrete summation '^2 i with the continuum limit 
^ J d9. We have checked by explicit computation^ that the term in curly brackets describing the correction to 
Gork'ov approximation is irrelevant within the approximation scheme adopted in previous Section and below, so we 
shall drop it from now on.[ This agrees with the fact that the holon scattering time behaves like ~ J 1-4 / 3 , a posteriori 
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a much longer time with respect to the electron scattering time, triggered by gauge fluctuations on spinons.] Then, 
the u integration can be performed exactly, after the trivial qq integration. The gii integration gives 



A " " I {X\\ — X°V F ) Xll — X V V F 



(54) 



where the last approximation is valid (in the weak sense) in the limit Av F x° 3> 1. Setting A(xq) = ^ v k p F x u J , the q± 
integration gives 



dq±e 



A(x ) 



A/A(x ) 



A/A(x ) 



dye lA{xo)x±v e^ v -A(x ) 



ji 1 A(i ) 



(=2 



(55) 



Collecting all pieces, the holon Green's function in the scaling limit and for Aa;o 3> 1 can be written as: 



A(x )fcF f ei^W 2 ^") 



A 



(16 



A{6) af3 e 



ikFXu (f?) 



7T7 7; Q(Xq) H T— „ 6( — X ) 

xi\{6)-x°v F y ' x\\ (9) + x°v F y ' 



e if*A^dx^ 



(56) 



C. Electron Propagator 

Next, we compute the electron Green's function using equation l|47|) . The Gor'kov terms in the ?p and the z 
correlators cancel against each other and the gauge field fluctuations act only on the gauge-invariant spinon correlator. 
Within the approximations used in previous Section, one can easily verify that this correlator in the scaling limit in 
x-space behaves as square root of the fi-propagator (|40(l . We perform now the Fourier transform of the electron 
propagator for momenta close to the Fermi surface, in a sector laballed by the angle 77, G a (uj, (it/2, tt/2) + n(i])kp +q), 
for small lo and q. We integrate over 9 using the following: 

Lemmaml Let f(9,x) be a smooth function, then in the large distance limit |a?| ^ A~ : we have 

J d9e lk -^- n ^ s f{9, x) ~ ^/(t?, x)S A -i (x A rt(r?)), (57) 

where 5\-i denotes an approximate <5-function of width A . 
Setting x = \x\n((f>) we approximate 

*A-i(£ATJ(77))~prp(0-»7)+*(0-»7 + 7r)]. (58) 
\x\ 

One can easily perform the (^-integration; the remaining integration over space-time variables is done as in previous 
Section, namely by saddle point approximation for |x| in the limit x° 3> \x\ and by principal part evaluation and scale 
renormalization for x°. The final result is 

2' 2J + fi ^ kp + 



s(v)z 



Jq-n{r])\x a (0) \ _ 



-iq-n(rt)\x c (0)\ 



V F 



d\x c \ 
d\xo\ 



(0)9|| 



v F 



d\x c \ 
d\xa\ 



(0)9|| 



where S(jj) is the angle dependent part of the wave function renormalization constant 



1 ~ 7^^ cos ^ + sin ( ? /)) 



(59) 



(60) 



This angle-dependent spectral weight is demonstrated in Fig. In (|59|l Z is the wave function renormalization 
constant averaged over the FS; writing the tomographic momentum cut-off as A = with A > 1 and taking into 
account the definition of Qq we obtain 



(61) 
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The rcnormalizcd electron self-energy E is given by 

/ . T 

E = !) s Wm; — ic— , (62) 

where we have reintroduced the spinon velocity v s previously set equal to 1. From H59JI we can immediately read off 
the inverse scattering time T for the electron: T = — ImE. 



D. Fermi Surface and Electron Resonance 



We now make the following 

Assumption FS: the neglected short - range attraction between spinon and holon renormalize the real part of E 
exactly to 0, so that the electron exhibits a Fermi surface. 

If Assumption FS holds, one might conjecture that this is due to a mechanism somewhat analogous to the one which 
renormalizes to zero the mass of the 2-fermion bound state ("pion") of massless QED3, whose constituent fermion are 
dynamically massive^ This cancellation between mass and self-energy attractions is there triggered by a symmetry 
principle, whose analogue in our scheme would require further investigation. One should remark that our treatment 
of the problem makes resemblance to the one discussed in Ref. Eil in the SU(2) x U(l) slave boson approach and in 
fact it yields a similar structure for the FS, although the scattering time is rather different. 

Under Assumption FS, finally, for q = and uj > small we find the structure (JSJ) with the replacement Z — * S(if)Z. 

This structure shows that the gauge fluctuations are able to bind together spinon and holon into a resonance for 
low energies and momenta close to the Fermi momenta, but with a wave function renormalization constant which 
depends both on the point of the FS, due to S(ri), and the temperature. In particular Z ~ T 1 / 6 , so Z vanishes 
if formally extrapolated to T = 0. This implies a peculiar non-Fermi liquid character for this system of "electron 
resonances" . However, a real extrapolation to T = cannot be done because the \x\ - saddle point is only dominant for 
T > x m sQo- The system therefore appears to fit naturally within the scheme of Unstable Fixed Points (UFP) outlined 
by Anderson^ There it is argued that in general in the Renormalization Group (RG) formalism, starting at high 
temperature and energy and integrating out high frequencies one derives a low temperature - low frequency model. 
However the system does not always flow smoothly under RG to T = 0. It might develop a tendency to approach at 
intermediate temperatures an infrared UFP. The composite holon-spinon system discussed above yielding the "electron 
resonance", might be a UFP in the temperature range of validity of our approximations. The angular dependence 
of the wave function renormalization S(rj) yields a reduction of the spectral weight outside the reduced Brillouin 
zone, in qualitative agreement with ARPES experiments in underdoped cuprates (see Ref. |2^ for a similar situation 
in the SU(2) x U(l) slave boson approach). In fact, the intensity measured in ARPES experiment, is proportional to 
ImG(w, k)n(uj); denoting by I(k) the integrated intensity along the "electron FS" we have a contribution to / due to 
the "electron resonance" given by 

I{n{6)k F ) ~ S(6)Z. 

The factor S{9) is peaked around 9 = ^ (for the FS near (?, ^)) and it is substantially reduced on the opposite 
side (see Fig. 01). 

One might try to extend the analysis performed above for momenta q around the holon FS in a shell of thickness A, 
by including the contributions of momenta q outside the shell, but still smaller than the U.V. cutoff of the continuum 
model. If in this contribution one tentatively neglects the effect of gauge fluctuations which give rise to incoherent 
component to the "electron" Green function, the essential features of our present consideration will remain; a detailed 
analysis is in progress. 



VI. COMPARISON WITH EXPERIMENTS 



A. MIC of In-plane Resistivity 



Let us now summarize the main results of Section IV, useful to derive a formula for the in-plane resistivity p, for 
comparison with the experimental data. 

We have shown that (under the stated approximations) the gauge fluctuations exhibit a typical scale, a sort of 
anomalous skin penetration depth, Qq 1 ~ (J- 2 /3y-i/3 _ j n ran g e f temperature identified by m 2 s > — > m s Qo 
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the gauge field couples a spinon-antispinon pair into a magnon resonance on a scale triggered by Q 
exhibits a complex mass term, M, of "relativistic" structure: 



The resonance 



M = 2, 



T 
ic—, 
X 



(63) 



where c ~ 3.3, whose imaginary part appears as a consequence of the dissipative nature of gauge fluctuations for 
energies smaller than T. 

The residue of the complex pole in the magnon resonance correlator is also T-dependent and behaves like 
Qo^Mk) 1 / 2 , where K is the Landau damping. Using j s ~ <90 ~ Qo^l and the Kubo formula for spinon conduc- 
tivity one obtains: 



1/2 



V5 



iarctanf^ 



(64) 



where £ ~ |^ ln<5|- 1/2 , A T ~ (x/Tc) 1 / 2 . 

For the holons one can borrow a computation performed diagrammatically in Ref. for a Fermi liquid interacting 
with a gauge field exhibiting Rcizcr singularity. Adding, via Matthiessen rule, the contribution of impurities one finds: 



Ph ~ SKepTimp)- 1 + (— ) 4/3 L (65) 

For small S, -j we have p s >> ph, so the spinon contribution dominates the physical resistivity in the Ioffe-Larkin 
rule. For low T, p s ~ -ip, thus exhibiting an insulating behaviour, for T > X m s one finds p s ~ T 1 / 4 , thus showing 
a metallic behaviour. From formula (|64|l a MIC is thus recovered decreasing the temperature, as shown in the 
experiments discussed in Sect. IC. This crossover is determined by the interplay between the AF correlation length £ 
and the thermal de Broglie wave length When At < £ the "peculiar" localization due to SRAFO is not felt and 

a metallic behavior is observed. In the opposite limite At > £ we find the insulating behaviour (but due to the gauge 
interaction , p s ^ e^) ; a behavior found for a "standard" localization). The doping dependence is rather weak, due 
to a delicate cancellation in the doping dependence of the dimensionless variable 

- T ) = ^?rws\' (66) 

which controls p s , sec (|64() . Our formula for p s has essentially no free parameters except for an overall resistivity 
scale. The only parameter 0(1) used in our numerical calculations is the coefficient r in the parametrization 

which one can fine-tune by using e.g. the minimum of resistivit y fo r some fixed doping. The entire set of curves 
p(S, T) arc then completely determined. As shown in Fig.l of Ref. I48L the agreement with experimental data is really 
good, and the MIC temperature goes down, as the doping increasing. If there were no logarithmic correction in our 
derived spinon mass, there would be no doping dependence AT ALL for the MIC temperature. 



B. Universal Normalized Resistivity 



Now we consider a more subtle prediction following from our theoretical treatment. 

As mentioned in Sect. IIIA, an inflection point T* has been observed in heavily underdoped cuprates at a higher 
temperature, where dp/dT has a maximum. Such an inflection point can be tentatively identified with the pseudogap 
temperature. Such an inflection point also appears in our derived in-plane resistivity formula (we still denote it by 
T*), and the "relativistic" structure of the mass term is responsible for it. We find that T* corresponds to y ~ 3.4 
Moreover, there is another inflection point in experimental data at low temperatures below which the resistivity 
exhibits an approximate logarithmic temperature dependence^ We propose to identify these two inflection points as 
the upper and lower bounds for the validity of our approximation and to approximately identify the 5 — T parameter 
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region corresponding to the "pseudogap phase" of cupratcs. Above the upper inflection point, the system enters the 
"strange metal" phase to which a separate paper is devotediS Below the lower inflection point, the system also crosses 
over to a new phase whose properties have to be explored. 

Neglecting the contribution, as justified at low T, from our formula (|64fl we notice that if we define the 
"normalised resistivity" 



p(T*) - p(T MIC ) 

this is a universal (i.e. doping independent) function of the variable y eq. I|tjtj|) . where Tju/c denotes the minimum of 
p s and one finds it corresponding to y ~ 1.7. 

This curve has been noticed in the YBCO^ data and quantitatively similar "universal curves" have been observed 
also for LSCOfS BSLCO, BSCOpii [In these last references a different definition of T* was used, based on deviation 
from linearity of p, not directly accessible to our approach, and therefore not permitting a direct comparison with our 
formula. A rough estimate, however, gives for that T* a value dependent on the material, but approximately twice 
our definition of T*]. Our formalism, on the other hand, explains in a neat way their universality character. In Fig.0] 
we plot the calculated normalized resistivity p n to be compared with the corresponding experimental curve on LSCO 
and YBCO that we extracted from the data of Takagi et al£ and Trappcniers et al.j2& see Fig. |3J We did not make 
any attempts to reconcile the calculated and observed location of MIC temperature which may depend on factors, 
not included in our consideration but the the universal character of the normalized resistivity is an explicit prediction 
of theory in agreement with experiment. 

Also, the recently experimentally observed a — b asymmetry in the conductivity of LSCO at low temperatures^ 2 ^ 
has a natural explanation in our framework as due to the anisotropy of the MIC temperature in a — b directions. A 
detailed explanation will be given in a separate communication^ 



C. Out-of-Plane Resistivity 

Let us summarize the results of Section V needed to compute p c . In the temperature range m 2 < — < m s Qo, 
the gauge fluctuations couple spinon and holon close to the FS into an "electron" resonance with scattering rate T 
proportional to the inverse life-time of the magnon, hence 




r = -Im Wm 2 -ic-~{ :,tsTi ] t (68) 



The wave function renormalization is the product of a term S(9), varying along the FS, where 9 is the angle labelling 
the direction from the center, inherited from the Dirac structure of the holon action, and a T-dependent term Z . This 
in turn is a product of a term proportional to the "magnon" renormalization constant Zq and a term coming from 
integration over fluctuations of holons, gaussian for those along and linear for those perpendicular to the FS, with a 

3 /2 

scale set by gauge fluctuations, hence giving a contribution ~ Q Q < 

Z ~ ZnQT ~ y/6m a Q . (69) 

To compute p c we average the angular dependence of S(9) and insert l|68|) and (|69J) in K-J's formula ea. lllOl) . 

It follows from (|68I1U|) that for low T, p c ~ T -1 and for higher temperature, if the first term in (|10|l still dominates, 
p c «~ T" 1 / 2 with a coefficient independent of t c . These features reproduce qualitatively the behavior observed 
experimentally in several materials (LSCO, YBCO... ) in the "pseudogap phase" including the rounded knee cited in 
Sect. IIIB, which corresponds to the above change of temperature dependence. 

As a consequence of K-J's approach p c at low T appears to give a direct test for the scattering rate of the "electron" 
in the pseudogap phase. The "metallic" contribution of the second term is important only at higher temperature, 
where it scales as T 1 / 6 , causing a further flattening of the p c (T) curves or possibly a minimum. Apart from an 
overall scale, having already fixed with p a b the variable x m s; our formula has only one free parameter, the scale of 
Z, i.e. essentially the scale A controlling the cutoff on momenta perpendicular to the FS, A = kp/X and weighting 
the "metallic" contribution. This parameter should be a somewhat large number and might be roughly estimated by 
fitting p c for one doping concentration. For other dopings the T-dependence behavior of p c is then derived and as one 
can see from Fig. the theoretical results are in good agreement with experimental data^ 

Having an explicit theoretical dependence on S and T for both p c and p a b one can further analyze the anisotropy 
ratio — . The derived temperature dependence of this ratio is shown in Fig. [7] this ratio clearly saturates at low T, 

Pab 
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since both p c and p scale as y but at higher temperature, in the "metallic" region for in-plane resistivity, the ratio 
decreases like T -1 / 4 . Again this behavior is qualitatively consistent with the experimental data in the "pseudogap 
phase"^* 4 ^ as shown in the inset of the same figure. 



D. Hidden MIC in Superconducting Cuprates and Magnetoresistance 

The techniques developed in previous sections are useful to compute other observables, like the transverse in-plane 
magnetoresistance and the 63 Cu spin-lattice relaxation rate. The calculation of magnetoresistance is outlined in Ref. 
I49I therefore here we only briefly review the results. 

The basic underlying hypothesis is that suppressing superconductivity by applying a magnetic field, in supercon- 
ducting underdoped samples one recovers the normal-state "pseudogap phase" . 

A magnetic field H perpendicular to the plane then modifies the gauge effective action in two ways: 1) Via a 
minimal coupling it induces a shift A — > A — eA em in the spinon term and A — > A + (1 — e)A em in the holon term, 
where e is the spinon effective charge and A em is the vector potential corresponding to the applied uniform static 
magnetic field H. In a mean-field treat ment t he effective charge should be chosen as to satisfy the Ioffc-Larkin rule 
for diamagnetic susceptibility (see Refs. I49l74() . Therefore e ~ Xh/x- 2) The presence of H induces a parity-breaking 
Chcrn-Simons term in the holon action (ah(H) /2ir)A° 'eijd 1 A? , where <Jh{H) is the holon Hall conductivity. Since Aq 
is short-ranged, with a gap 7 = v + uj p (see ea. (|22ll ). it can be integrated out first yielding an effective renormalization 

of the diamagnetic susceptibility in the transverse action: x ~ * x{H) = x + as discussed in Ref. |3 This effect 

is however subleading at low T. 

Under the approximations of Sect. IV the result of effects 1) and 2) can be summarized by a modification of the 
"relativistic" mass term of spinon: 



"-"^-ir-'Gssrw)- (70) 

[A technical comment: the minimal coupling in the FSF path-representation of G(x,0\F) produces a term 
exp[ie J Q dXX Jq ds' ds"{(j) l {s') — 2p l ){(f>> \s") — 2p>)eijH] , see ea. (|2"tj|) . and evaluating the <f> integral in gaussian 
approximation this term yields a contribution e 18 to the Q correlator in ea. (|31[l . This finally is responsible 

for the shift of the last term in the square root in eq. I|70jl] The limits of validity of the \x\ saddle point become 
x(H)Qo\M(H)\ < T, lm(M(H)) < m 2 s , which for the range of physical parameters considered here {H < 100 Tesla) 
gives a temperature range still lying between a few tens and a few hundreds degrees. Hence, to conclude, the presence 
of H modifies p s via the cyclotron effect, by reducing the damping from to ^fm — ^[qL • 

This reduction makes the thermal de Broglic wavelength At longer, so the MIC occurs at a higher temperature 
w.r.t. the system at H = 0. The external magnetic field then reveals the MIC originally hidden in the superconducting 
samples. 

Furthermore the shift of the minimum of p causes a strong positive transverse magnetoresistance (MR) at low T, 
as in fact experimentally seenj2Si21 an effect missing in previous theoretical treatments^ At higher temperatures, in 
the region where dissipation dominates, the shift of diamagnetic susceptibility due to the Chern-Simons term induces 
a reduction of resistivity, a tendency contrasted by the classical cyclotron effect on holons, taken into account in the 
Boltzmann equation approximation. One then has two possible types of MR curves: one is always positive but it 
exhibits a knee below the crossover temperature between the mass gap and the dissipation dominated regions (See 
Fig. 3 in Ref. l49j) . This behavior can be compared with the one observed in LSCO reported in Ref. an d one 
finds a reasonably good agreement. If, on the contrast, the quantum effects related to (Jh{H) are sufficiently strong, 
a minimum develops, eventually leading to a negative MR in some region around it. The MR scales quadratically 
with H (See Fig. 2 in Ref. ) in agreement, in particular, with data on LSCOj^l away from the doping 6 = 1/8 
where the stripe effects dominate. [In the explicit formula one should also take into account the modification induced 
by H in the contribution of the Landau damping to Zf \J~k ~^( x (jj") QcT 3 ~ cIH^Qq ) x / 2 where d is a new constant 
- /''(Ce^/4) in fact 

roughly estimated, together with a parameter coefficient of H in cr(iJ), by comparison with an 
experimental curve at some doping.] 

Finally we notice that in Zn-doped superconducting samples of BSLCO the MIC become observable upon increase 
of Zn doping (when a magnetic field suppresses superconductivity) and it shifts to higher temperature as the level 
of Zn-doping increases^ This effect is qualitatively consistent with our picture. In fact, the Zn-doping disturbs the 
AF background, so making the AF correlation length shorter, therefore shifting the MIC temperature up, although 
we are not able, at the moment, to make a quantitative estimate of this shift. 
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E. Spin-Lattice Relaxation Rate 

We turn now to the spin-lattice relaxation rate for the Cu-sites, (tjt-) 63 , see Ref. 0. This can be theoretically 
computed using the Kubo formula: 

^ = lim / Al^l 2 ^^, (71) 



T\T J oj 

where A(q) is the hyperfme field and Xs(q^) the spin susceptibility. For the Cu sites the hypcrfinc field A(q) is 
peaked around Qaf = (tt, n), thus probing the AF spin fluctuations. The electron spin field S(x) = ^c{x) is related 
to the spinon and holon fields by 

S(x) - (f - H*H(x))e^ AF - s n(x). 
Approximating H* H by its mean field value 5 and using the Lehmann representation one finds, for small q: 

lim lm Xs(QAF+q,Uj) _ Im f°° ^ f ^ (1 _ §)2 < fa, ^ (Q) > ^.s_ (?2) 
U->0 LU Jo J 

The |x| and x° integration are performed as in Sect. IVC. Assuming a cutoff for the \q\ integration in l|7fl ) given 
by the inverse anomalous skin depth, Qq, and using the smoothness of A(q) at this scale we derive 

d6 f d\q\\q\ \A(q)\ 2 e l ^ s ^ cose ~Q 2 J {Ce™/ 4 ). (73) 

Numerically one finds ReJ (ce 47r/4 ) = a and Im J (ce i7T / 4 ) = b with a/b ~ 0.1. Plugging J73J) and gUJ) in J72J) one 
obtains from the Kubo formula l|7f ll 

(T lT r - (1 - SrVS\MrHacos(^) + 6sin(^)) (74) 

For low T 7 Tfkp ~ a + bT and for higher T one finds ~ T" 1 / 4 ; therefore the spin lattice relaxation rate ( jr^) 
on Cu- sites cxnibits a maximum and an inflection point at higher temperature, as observed in YBCO underdoped 
samples 

If a would be 0, then would be proportional to the spinon conductivity <r s , and the maximum and the inflection 
point would be at the same temperature of the MIC and T*, respectively. However, due to the a term in ifT^ they are 
shifted. In particular, the inflection point is found at a lower temperature Tq, in qualitative agreement with the fact 
that experimentally the pseudogap temperature deduced from spin-lattice relaxation rate is lower than that derived 
from the resistivity measurements S 

We end this Section by remarking that preliminary calculations are giving also encouraging results, when compared 
with the experimental data, for the electronic AC conductivity at small and the electronic specific heati^ 

VII. CONCLUDING REMARKS 

To summarize we have presented in this paper the calculation of physical quantities like the in-plane and out-of- 
plane resistivities, spin-lattice relaxation rate, etc within the spin-charge gauge field approach, and compared the 
theoretical results with experimental data in the pseudogap phase with a very good agreement. In particular, we have 
elucidated the origin of the MIC in the non-superconducting cuprates as well as the MIC in superconducting samples 
when a strong magnetic field suppresses the superconductivity. In our view, this striking phenomenon is an intrinsic 
property of the pseudogap phase which can shed light on other puzzles in this regime. We are still in the process of 
studying the diversified properties using our approach in this interesting phase. 

There is some skeptism w.r.t. the gauge field approach in general, mainly because of the strong interactions among 
the constituent particles. Our attempt in this direction, at least, gives some encouraging signal: if the underlying 
physics is grasped by the treatment, and appropriate non-perturbative tools are employed, there is a fair chance to 



21 



correctly describe the puzzling phenomena in the strongly-correlated systems. The treatment is not "rigorous" in the 
mathematical-physics sense, but still acceptable by "theoretical-physics" standard. Needless to say, the final word 
belongs to experiments, verifying all consequences of the theoretical interpretation. 

The gauge field approach provides us with a non-trivial picture in strongly correlated two-dimensional systems. 
Unlike the one-dimensional systems where the spin and charge are fully separated, and three-dimensional systems 
where the spin and charge are confined, the spin and charge in two-dimensional systems appear "separated" in their 
scattering against gauge fluctuations, while being bound into "electron" resonance at low energy- momentum scale. In 
particular, in the "pscudogap phase" the presence of 7r flux and AF Neel background makes the system " relativistic" 
with linear dispersion for the holons and spinons " massive" due to interaction with vortices attached to slowly moving 
renormalized holons. If the system were truly "relativistic", we would have spinon-holon confinement. However in 
the actual system the presence of a finite Fermi surface breaks the "bootstrap" symmetry and gives rise to the very 
peculiar Reizer singularity, producing the binding force between spinon-antispinon and spinon-holon. The physical 
consequences of this non-trivial picture have to be further explored. 
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FIG. 1: Feynman-Schwinger-Fradkin representation 



FIG. 2: Tomographic decomposition of the Fermi Surface patch with square boxes of size A. 
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FIG. 3: Angle-dependent spectral weight of the electron propagator. The thick lines close to (±7r/2, ±7r/2) represent the region 
of FS with spectral weight larger than 1/2 for 8 ~ 0.05 
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FIG. 5: Temperature dependence for p n in underdoped LSCO ( Extracted from Ref. 0) and YBCO (Extracted from Ref. l3^) 
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FIG. 6: Calculated temperature dependence of the out-of-plane resistivity (in arbitrary units) for different doping concentra- 
tions: 5=0.05 (full line), 5=0.07 (dashed) and (5=0.09 (dotted). Inset shows experimental data on LSCO, extracted from Ref. 
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FIG. 7: Calculated temperature dependence of the resistivity anisotropy ratio as a function of temperature for different doping 
concentration: (5=0.05 (full line), S =0.07 (dashed) and 5=0.09 (dotted). Inset shows corresponding experimental data on 
LSCO, extracted from Ref. [fit] 



